AD-A169  <78 
UNCLASSIFIED 


IDENTIFICATION  OF  DAMP I NO  PARAHETERS  OF  A  SOUEEZE-FILH 
BEARINQ(U)  SUSSEX  UNIV  BRIGHTON  (ENGLAND)  J  8  ROBERTS 
NAV  86  DAJA43-84-C-8844 

F/G  28/11 


1/1 


1.1 


.6 


AD-A169  678 


lr  v  ^  w '  r1  wj  '.Rl'R  1 
yNcLAss^fcifeu 


SECURITY  CLASSIFICATION  OF  THU  PACE  (Whom  Dm m  Entered) 


■  Tt to  *y  '3  ’T M  ^  FJM  jr  JVjyj  A  |R.  r.'  rj  rvr.'r: 

R&D  4005A-AN 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  REPORT  NUMBER 


[2.  GOVT  ACCESSION  NO. 


>.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (end  Subtitle) 

Identification  of  Damping  Parameters  of  a 
Squeeze-Film  Bearing 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Technical  Report 
Oct  84  -  Oct  85 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR**; 

Dr.  J.B.  Roberts 


CONTRACT  OR  GRANT  NUMBER**; 

DAJA45-84-C-0044 


».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Sussex 
U.K. 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

61102A 

1L161102BH57-06 


11.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

USARDSG-UK 

Box  65,  FPO  NY  09510-1500 


12.  REPORT  DATE 

May  1986 


IS.  NUMBER  OF  PAGES 

46 


U.  MONITORING  AGENCY  NAME  4  ADDRESS*!/  dl  He  rent  from  Controlling  Oltlce) 


IS.  SECURITY  CLASS,  (ol  thle  report.) 

Unclassified 


ISa.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  ill la  Report) 


Approved  for  Public  Release;  distribution  unlimited 


19.  KEY  WORDS  (Centlnuo  on  rerera*  alda  II  nacaaaarr  and  Idontlty  by  block  number) 

Vibration  Control;  Squeeze-film  bearings;  Damping  measurements; 
Parametric  identification;  Inertial  effects 


O- 

o 

IjU 
_ I 

uz 


20k  ABSTRACT  fGrtltaa  an  rarer**  *M>  R  nmemaom ry  mad.  Ideality  by  block  nuaibar; 

This  report  describes  the  results  obtained  from  a  combined  experimental  and 
theoretical  programme  concerned  with  the  determination  of  the  damping  and 
inertial  coefficients  of  a  squeeze-film  bearing,  and  the  dependence  of  these 
coefficients  on  eccenticity  ratio  and  frequency  of  vibration. 

The  experimental  results  were  obtained  from  a  new  experimental  rig,  in  which 
the  squeeze-film  journal  is  mounted  on  a  beam  of  continuously  variable  length. 


■  u  hi  ■  •  •Al-  ” " ”  ~ ^  -  v  J.v.’  tv  ■  -■  .■ * •'.■  tv.-.v .y 

SECURITY  CLASSIFICATION  OF  THI»  PAOtfWhan  Data  lm«w<) 

20.  Oontd. 

By  means  of  a  transient  test  technique ,  combined  with  signal  processing, 
based  on  a  parametric  identification  method,  estimates  of  the  damping  and 
inertial  coefficients  have  been  obtained  for  static  eccentricity  ratios  in 
the  range  0.0  to  0.5,  and  for  frequencies  of  vibration  ranging  from  50  Hz 
to  180  Hz.  These  coefficients  are,  in  general,  much  higher  in  value  than 
the  conventional  short-bearing  theory  for  squeeze-film  bearings  predicts 
(by  a  factor  of  10,  in  some  cases)  suggesting  serious  inadequacies  in  the 
available  theoretical  prediction  procedures. 

In  an  attempt  to  explore  possible  theoretical  deficiencies,  with  respect  to 
the  evaluation  of  inertial  coefficients,  numerical  theoretical  predictions 
of  these  coefficients  were  obtained  for  finite-length  bearings,  with  various 
L/D  ratios  (where  L  is  the  length  of  one  bearing  land,  and  D  is  the  diameter 
of  the  bearing) .  Hie  technique  used,  based  on  a  finite  difference  solution 
of  a  two-dimensional  partial  differential  equation  for  the  "inertial 
pressure"  distribution,  was  found  to  give  good  agreement  with  some  other 
theoretical  results,  obtained  by  much  more  complex  methods.  However,  the 
effect  of  the  finite  width  of  the  bearing  was  found  to  be  weak,  for  the 
squeeze-film  geometry  of  interest  here  -  i .  e . ,  the  short-bearing  theory  was 
found  to  be  an  adequate  approximation.  Hius,  this  theoretical  investigation 
has  not  revealed  the  basic  source  of  discrepancy  between  the  measured 
coefficients  and  the  short-bearing  theory,  based  on  conventional  assumptions 
regarding  the  fluid  flow. 

This  report  describes  the  experimental  and  theoretical  aspects  of  the  present 
investigation  in  two  parts,  referred  to  as  part  one  and  two,  respectively, 
followed  by  some  concluding  remarks  on  the  investigation  as  a  whole. 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWhan  Data  Bnltnd) 


IDENTIFICATION  OF  DAMPING 
PARAMETERS  OF  A  SQUEEZE-FILM  BEARING 


Final  Technical  Report 


by 


Dr.  J.B.  Roberts 


May  1986 


United  States  Army 

EUROPEAN  RESEARCH  OFFICE  OF  THE  U.S.  ARMY 
London,  England. 


CONTRACT  NUMBER  DAJA45-84-C-0044 


University  of  Sussex, 


Approved  for  Public  Release;  distribution  unlimited 


Summary 

This  report  describes  the  results  obtained  from  a  combined 
experimental  and  theoretical  programme  concerned  with  the 
determination  of  the  damping  and  inertial  coefficients  of  a 
squeeze-film  bearing,  and  the  dependence  of  these  coefficients 
on  eccentricity  ratio  and  frequency  of  vibration.  i  , 

The  experimental  results  were  obtained  from  a  new 
experimental  rig,'  in  which  the  squeeze-film  journal  is  mounted 
on  a  beam  of  continuously  variable  length.  By  means  of  a 
transient  test  technique,  combined  with  signal  processing^, 
based  on  a  parametric  identification  n^ethod,  estimates  of  the 
damping  and  inertial  coefficeints  have  been  obtained  for 
static  eccentricity  ratios  in  the  range  0.0  to  0.5,  and  for 
frequencies  of  vibration  ranging  from  50  Hz  to  180  Hz.  These 
coefficients  are,  in  general,  much  higher  in  value  than  "the 
conventional  short-bearing  theory  for  squeeze-film  bearings 
predicts  (by  a  factor  of  10,  in  some  cases)  suggesting  serious 
inadequacies  in  the  available  theoretical  prediction 
procedures . 

In  an  attempt  to  explore  possible  theoretical 
deficiencies,  with  respect  to  the  evaluation  of  inertial 
coefficients,  numerical  theoretical  predictions  of  these 
coefficients  were  obtained  for  finite-length  bearings,  with 
various  L/D  ratios  (where  L  is  the  length  o^E  one  bearing  land, 
and  D  is  the  diameter  of  the  bearingjv  ThW  technique  used, 
based  on  a  finite  difference  solution  of  a  two-dimensional 
partial  differential  equation  for  the  "inertial  pressure" 
distribution,  was  found  to  give  good  agreement  with  some  other 
theoretical  results*  obtained  by  much  more  complex  methods. 
However,  the  effect  of  the  finite  width  of  the  bearing  was 
found  to  be  weak,  for  the  squeeze-film  geometry  of  Interest 
here  -  i.e.,  the  short-bearing  theory  was  found  to  be  an 
adequate  approximation.  Thus,  this  theoretical  Investigation 
has  not  revealed  the  basic  source  of  discrepancy  between  the 
measured  coefficients  and  the  short-bearing  theory,  based  on 
conventional  assumptions  regarding  the  fluid  flow. 

This  report  describes  the  experimental  and  theoretical 
aspects  of  the  present  investigation  in  two  parts,  referred  to 
as  part  one  and  two,  respectively,  followed  by  some  concluding 
remarks  on  the  Investigation  as  a  whole. 
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NOTATION 


true  non-dimensional  start  displacement  for 
free-decay  data 

direct  fluid  damping  coefficients,  for  the  r 
and  s  directions,  respectively 

non-dimensional  direct  fluid  film  damping 
coefficients,  for  the  r  and  s  directions, 
respectively 

radial  clearance  between  journal  and  bearing 

direct  fluid  film  inertial  coefficients,  for 
the  r  and  s  directions,  respectively. 

initial  displacement 

non-dimensional  initial  displacement  (*  d/c) 

mean  square  error  between  fitted  theoretical 
decay  curve  and  experimental  data 

beam  stiffness,  in  the  r  and  s  directions, 
respectively 

land  length  of  bearing 

effective  first  mode  mechanical  mass 

hydrodynamic  mass  (crr  in  radial  direction, 
css  in  transverse  direction) 

journal  displacement  in  the  radial  direction 
(i.e  in  the  direction  of  the  attitude  line) 

radius  of  the  journal 

Reynolds  number  (  =  puc2/n) 

journal  displacement  in  the  transverse  direction 
(i.e.  in  a  direction  perpendicular  to  the 
attitude  line) 

t  ime 

guessed  start  time 

true  start  time 

non-dimensional  displacement 

measured,  free-decay  record 

analytical  solution  to  equation  (4) 

objective  function,  defined  by  equation  (11) 


non-dimensional  damping  parameter 

value  of  y*  ,  in  the  absence  of  the  fluid  film 

non-dimensional  stiffness  parameter 

value  of  S*,  in  the  absence  of  the  fluid  film 

eccentricity  ratio 

static  eccentricity  ratio 

parameter  defined  by  equation  (10) 

absolute  viscosity  of  the  squeeze-film  fluid 

density  of  the  squeeze-film  fluid 

sampling  interval 

difference  between  tg  and  t? 

frequency  of  vibration 

natural  frequency  of  undamped  vibration 
(u  =  uro,  for  radial  motion,  and 
u>  =  uso ,  for  transverse  motion) 

non-dimensional  time 


SUMMARY 


A  technique  for  obtaining  estimates  of  the  direct  damping 
and  inertial  coefficients  of  a  squeeze-film  bearing  is 
described.  This  involves  applying  parametric  identification 
and  optimisation  techniques  to  digitised,  free-decay 
experimental  displacement  records.  The  experimentally  obtained 
coefficients,  derived  by  this  technique,  were  found  to  be 
significantly  higher  in  magnitude  (in  some  cases  by  a  factor 
of  about  10)  than  the  corresponding  values  derived  from 
conventional  short-bearing  theory,  and  to  be  virtually 
independent  of  the  frequency  of  vibration. 

1 .  INTRODUCTION 

The  squeeze-film  damper  has  proved  a  useful  means  of 
eliminating  instabilities  and  reducing  vibration  in 
rotor-bearing  systems,  such  as  occur,  for  example,  in 
aero-engines  (e.g.  see  Refs  1  and  2). 

For  squeeze-film  dampers  with  a  centralising  spring  it  is 
possible  to  model  their  dynamic  behaviour  in  terms  of  damping 
coefficients,  which  can  be  estimated  theoretically  from  a 
solution  of  Reynolds  equation  (e.g.  see  Ref. 3).  Recently  it 
has  been  demonstrated  theoretically,  by  a  number  of  workers 
[4-10],  that  inertial  effects  within  the  lubricant  can  also  be 
significant  and  that  the  normal  damping  coefficient 
representation  should  be  augmented  to  include  acceleration,  or 
inertial  coefficients.  In  many  applications  these  latter 
coefficients  are  important  since  the  gap  Reynolds  number 
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Despite  these  advances  in  theoretical  understanding  there 
has  been  a  distinct  lack  of  supporting  experimental 
measurements,  particularly  with  regard  to  the  effect  of  fluid 
inertia,  and  agreement  between  theory  and  available 
experimental  results  is  generally  rather  inconclusive  (e.g. 
see  Refs.  11-14).  In  an  effort  to  improve  this  situation  a 
recent  experimental  investigation  was  carried  out,  using 
parametric  identification  techniques,  for  the  purposes  of 
obtaining  estimates  of  damping  and  inertial  coefficients  from 
transient  test  data  [15].  The  results,  which  covered  the 
Reynolds  number  range  0.5  <  Re  <  1.5,  indicated  that  the 
inertial  effect  was  significantly  greater  than  that  predicted 
by  available  theories. 

The  objective  of  the  present  paper  is  to  present  some  new 
experimental  measurements  of  damping  and  inertial 
coefficients.  The  same  basic  test  procedure  as  described  in 
reference  [15]  was  employed  but  the  data  was  obtained  from  an 
entirely  new  rig,  which  allowed  a  continuous  variation  of  the 


frequency  of  vibration.  Moreover,  the  squeeze-film  geometry 
used  in  the  present  investigation  was  closer  to  that  used  in 
typical  aero-engine  applications.  The  results  relate  to  a 
higher  frequency  range  (50Hz  to  180Hz)  than  that  in  [15];  the 
corresponding  Reynolds  number  range  is  1.05  <  Re  <  3.64. 

2  THE  EQUATIONS  OF  MOTION 

For  a  squeeze-film  damper  with  a  centralising  spring, 
small  motion  of  the  journal,  with  respect  to  a  fixed  bearing, 
can  be  described  by  a  pair  of  second  order  equations  of 
motion.  These  equations  are  conveniently  written  in  terms  of 
the  radial,  r,  and  transverse,  s,  displacements,  where  these 
are  measured  from  the  static  eccentricity  position.  For  free 
motion,  in  the  radial  direction  only,  the  equation  of  motion 
for  r  is  simply 

(m+crr ) f  +  (bg  +  brr ) t  +  krr  =  0  (1) 

where  crr  is  a  radial  inertial  {or  acceleration)  coefficient 
and  brr  is  a  radial  damping  coefficient  for  the  squeeze-film, 
bg  denotes  structural  damping  and  kr  is  the  stiffness  of  the 
centralising  spring;  m  is  the  mass  of  the  squeeze-film 
journal.  For  free  motion  in  the  transverse,  s,  direction  an 
equation  of  motion  similar  to  equation  (1)  applies;  the 
relevant  squeeze-film  coefficients  are  now  css ,  the  transverse 
inertial  coefficient  and  bSs  *  the  transverse  damping 
coefficient. 

Both  these  equations  of  motion  for  free  vibration  can  be 
written  in  a  standard  form.  As  a  first  step  it  is  convenient 
to  non-dimensionalise  by  introducing  the  non-dimensional  time 
variable 


t  =  uQt  (2) 

where  uQ  is  the  frequency  of  undamped  vibration  (uD  =  wor ,  for 
radial  motion),  and  uG  =  uos,  for  transverse  motion,  and  the 
non-dimensional  displacement 


y  i  *■  \  _  displacement  (at  t) _ 

'  initial  displacement  (at  t=0 )  *  ' 

The  standard  form  is  then 

X  +  y*  X  +  S*  X  =  0  (4) 

where  y*  is  a  non-dimensional  damping  coefficient  and  6*  is  a 
non-dimensional  stiffness  coefficient. 

Through  a  comparison  of  equations  (1)  and  (4)  one  readily 
finds  that,  for  the  radial  motion, 

brr  =  m  uro  £  ^pjr  -  yQ*j  (5) 


and 


err  ■  m  f  4*  “  1  1 


(6) 


where  yQ*  is  the  value  of  y*  in  the  case  of  no  fluid  film  in 
the  squeeze-film  damper  -  i.e.  the  structural  damping 
component  of  y* .  Similar  equations  for  bss  and  css  can  be 
written. 

Equations  (5)  and  (6)  provide  a  simple  means  of 
determining  brr  and  crr  (and  also,  of  course,  bss  and  css) 
from  free-decay  tests.  Thus,  for  example,  from  a  measured 
free-decay  in  the  radial  direction  only,  values  of  y*  and  S* 
may  be  found  by  fitting  equation  (4),  in  some  sense,  to  the 
experimental  data.  Similarly,  yQ*  and  uro  may  be  determined 
from  a  free-decay  test  without  lubricant  in  the  damper. 
Equations  (5)  and  (6)  can  then  be  applied  to  yield  estimates 
of  brr  and  crr . 

In  the  tests  to  be  described  later  the  centralising  spring 
was  realised  by  mounting  the  squeeze-film  journal  centrally  on 
a  beam,  built-in  at  both  ends.  Free-decay  motion  was  obtained 
by  applying  a  force  to  the  centre  of  the  journal,  to  give  the 
journal  an  initial  displacement,  and  then  releasing  the  force. 
For  such  a  system  many  modes  of  vibration  can,  in  principle, 
exist  simultaneously;  however,  it  is  easy  to  show  that  the 
first  mode  will  be  entirely  dominant  and  thus  that  an  equation 
of  motion  of  the  form  of  equation  (4)  is  appropriate.  In 
these  circumstances  the  mass  m  appearing  in  equations  (5)  and 
(6)  can  be  interpreted  as  the  effective,  first  mode  mass.  It 
can  be  determined  very  easily  by  a  direct  measurement  of  the 
beam  stiffness,  k,  and  the  undamped  natural  frequency  of 
vibration,  u0. 

For  the  case  of  a  short  full-film  squeeze-film  damper, 
which  is  of  Interest  here,  theoretical  expressions  of  b^,.  and 
bss  have  been  known  for  some  time  [16].  Corresponding 
theoretical  results  for  crr  and  css  have  been  given  by  Smith 
[17].  These  results  are  collected  together  in  Table  1.  It 
has  been  shown  recently  by  the  authors  [18]  that  Smith's 
theoretical  results  are  in  very  close  agreement  with  the 
results  of  the  more  elaborate  analyses  given  by  Reinhardt  and 
Lund  [8]  and  Szeri  et  al  [9] 

It  is  noted  that  these  theoretical  results  are  strictly 
valid  for  small  values  of  Re  (Re  <<  1),  whereas  in  the  present 
experiments  Re  >  1 .  However,  theoretical  investigations  into 
the  influence  of  finite  Reynolds  number  on  the  force 
coefficients  (e.g.  see  Refs  19  and  20)  suggest  that,  for  the 
squeeze-film  geometry  of  concern  here,  this  influence  is 
likely  to  be  small. 

3.  EXPERIMENTAL  FACILITY  AND  PROCEDURE 
3 . 1  Experimental  rig 


The  experimental  rig  is  shown  by  a  photograph  in  Fig. 1(a) 
and  as  a  line  diagram  in  Fig.  1(b).  The  journal  of  the 
squeeze-film  damper  was  non-rotating  but  free  to  vibrate,  as  a 
result  of  external  disturbances,  in  any  direction  orthogonal 
to  its  axis.  It  was  contained  within  a  circular  bearing  and 
supported  by  a  beam  of  circular  cross-section.  This  beam, 
which  provided  a  stiffness  in  parallel  with  the  squeeze-film 
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(i.e.  acted  as  a  centralising  spring)  was  clamped  rigidly, 
through  a  collet  connection,  to  support  brackets,  positioned 
symmetrically  on  each  side  of  the  journal.  The  effective 
length  of  the  support  beam  could  be  varied,  fairly  easily,  by 
varying  the  position  of  the  support  clamps.  In  this  way  it 
was  possible  to  vary  the  effective  spring  stiffness  of  the 
beam,  and  hence  the  natural  frequency  of  the  first  mode  of 
vibration,  over  a  fairly  wide  range.  In  addition  to  this 
adjustment,  screw  actuated  positioning  devices  allowed  the 
bearing  to  be  moved  independently  in  both  the  horizontal  and 
vertical  directions  (perpendicular  to  the  beam  axis) ;  thus  any 
desired  static  eccentricity  position,  of  the  journal  relative 
to  the  bearing,  could  be  readily  achieved.  Accurate  static 
positioning  was  obtained  by  outputting  displacement  transducer 
signals  to  the  computer  display  screen  whilst  the  operating 
screws  were  adjusted. 

The  bearing  consisted  of  two  plain  lands  separated  by  a 
central  circumferential  groove.  Lubricant  (in  this  case 
Tellus  RIO)  was  supplied  to  the  groove  by  a  pump  through  three 
2BA ,  3/4"  deep  holes,  spaced  120°  apart  around  the  bearing. 
No  end  seals  were  employed,  so  the  lubricant  simply  discharged 
freely  into  a  reservoir  before  being  recirculated.  The 
relevant  bearing  dimensions  were  as  follows: 


Bearing  land  length  9.0mm 

Journal  radius  68.0mm 

Radial  clearance  0.22mm 

Central  groove  width  4.0mm 

Central  groove  depth  0.5mm 


3 . 2  Experimental  procedure  and  data  acquisition 

In  the  experimental  work,  tests  were  carried  out  for 
natural  frequencies  of  vibration  of  50Hz,  80Hz,  128Hz,  150Hz 

and  180Hz  (it  is  noted  that  the  differences  between  and 

uso  were  found  to  be  very  small  in  practice;  this  distinction 
will,  therefore,  not  be  made  in  future  discussion).  In  each 
case  the  effective  spring  stiffness  of  the  beam  was  determined 
experimentally  by  simply  pulling  the  journal  horizontally, 
with  a  known  static  force,  and  measuring  the  resulting  journal 
displacement.  From  this  information,  and  the  measured  natural 
frequencies,  it  was  possible  to  determine  the  effective 
journal  mass,  m,  at  each  natural  frequency. 

For  each  of  the  natural  frequencies  chosen,  a  series  of 
test  was  conducted,  at  various  values  of  static  eccentricity 
ratio.  In  each  test  the  journal  was  pulled  through  an 
"initial  displacement"  to  an  "initial  position"  within  the 
clearance  circle,  by  a  length  of  wire  looped  over  a  pin 
attached  to  the  core  of  a  solenoid.  The  journal  was  then 
released  by  actuating  the  solenoid;  this  had  the  effect  of 
rapidly  retracting  the  pin  and  hence  freeing  the  wire.  The 
subsequent  transient,  decaying  oscillatory  motion  of  the 
journal  was  measured  by  proximity  capacitive  probes  clamped 
adjacent  to  datum  surfaces  on  the  vibrating  journal  and 
recorded  by  computer. 


In  the  present  series  of  tests  the  static  eccentricity 


position  was  always  on  the  vertical  radial  line  and  radial  and 
transverse  free  decays  were  measured  separately  by  pulling  the 
journal  in  the  vertical  and  horizontal  directions, 
respectively.  For  each  static  eccentricity  value  results  were 
obtained  for  (a)  the  journal  centre  the  required  distance 
above  the  bearing  centre  and  (b)  the  journal  centre  the  same 
distance  below  the  bearing  centre,  and  the  results  finally 
averaged . 

For  each  transient  test  it  was  possible  to  measure, 
simultaneously,  both  horizontal  and  vertical  motions.  It  was 
observed  experimentally  that,  as  one  would  expect  from  the 
theory,  a  perturbation  of  the  journal  in  one  direction  (either 
radial  or  transverse)  did  not  result  in  any  significant  motion 
in  the  perpendicular  direction.  In  fact,  motion  in  the 
perpendicular  direction  was  found  to  be  less  than  one  per  cent 
of  the  amplitude  of  the  motion  in  the  direction  of  initial 
displacement.  Thus,  in  the  main  body  of  the  experimental 
work,  only  the  motion  in  the  direction  of  the  initial 
displacement  was  measured  and  recorded. 

Throughout  the  tests  the  supply  pressure  of  the  lubricant 
was  maintained  at  14  kN/m2 ,  to  ensure  that  full-film 
lubrication  conditions  pertained  in  the  squeeze-film.  The 
temperature  of  the  lubricant  was  measured  at  the  outlet  of  the 
bearing  and  recorded  for  each  test.  The  corresponding 
lubricant  viscosity  was  then  obtained  from  a  calibration 
curve.  The  viscosity  changed  significantly  with  temperature. 
Thus  at  20°c,  n  «  0.0188  Ns/m2,  whereas  at  30°C,  n  «=  0.0122 
Ns/m2 . 

To  avoid  the  necessity  of  processing  and  graph  plotting 
the  decay  data  manually,  a  microprocessor  based,  purpose-built 
data  acquisition  system  was  developed  for  use  in  this  work; 
this  could  also  be  used  for  initiating  a  test.  An  initiation 
pulse,  generated  by  the  microprocessor  was  used  to  activate  a 
relay  and  hence  energise  the  solenoid  and  release  the  journal; 
this  pulse  also  simultaneously  started  analogue  to  digital 
conversion  (ADC)  of  the  journal  displacement  signal  coming 
from  the  rig.  The  sampling  interval,  At,  in  the  ADC 

conversion  could  be  varied,  under  software  control.  The 
equi-spaced  samples  of  the  free-decay  were  stored,  and  then 
transmitted  to  a  main-frame  computer,  for  subsequent 
processing. 

4  ESTIMATION  OF  PARAMETERS 

As  explained  in  the  preceding  section,  the  proposed 
procedure  requires  the  estimation  of  y*  and  s* ,  in  equation 
(4),  from  normalised,  experimental  free-decay  data.  The  data 
was  in  the  form  of  n  equi-spaced  digital  samples,  Xj , 

(i  =  0 , 1 , . . . , n-1 ) .  The  estimation  technique  employed  here  is  a 
combination  of  a  parametric  identification  method  (previously 
used  in  ref  [15])  and  an  optimisation  method. 

4 . 1  Parametric  identification 

In  ref  [15]  an  iterative  algorithm  for  identifying  the 
parameters  y*  and  S*  in  equation  (4)  was  developed,  based  on 
the  sequential  least-square  method  of  Detchmendy  and  Stridhar 
[21].  This  method,  also  used  in  the  present  work,  generates 


sequential  estimates  of  the  state  variables  X(t)  and  X(t)  and 
also  estimates  of  y*  and  S*  ,  at  the  times  iAt.  These 
estimates  "track"  the  experimental  data  and  generally  converge 
to  stable  values  as  time  increases.  By  applying  the  method 
iteratively,  improved  estimates  of  y *  and  6*  can  be  found. 
Thus,  in  the  first  iteration,  y*  and  S*  are  guessed,  and  these 
values  are  used  to  start  the  sequential  estimation 
computation.  The  estimates  obtained  at  the  end  of  the  data 
sequence  are  then  used,  in  place  of  the  initial  guess,  to 
start  the  next  iteration.  By  repeating  the  iteration  a  number 
of  times,  the  estimates  of  y*  and  S*  ,  at  the  end  of  the  data 
sequence,  are  found  to  converge  to  constant  values. 

4 . 2  Optimisation  method 

A  difficulty  with  the  method  just  described  is  that  the 
exact  time  of  release,  i.e.  the  true  start  time  of  the  free 
decay,  is  not  known  accurately.  This  is  because  there  is  a 
finite  time  required  to  energise  the  solenoid  and  retract  its 
pin.  This  delay  time  was  found  to  vary  somewhat  from  test  to 
test,  due  to  the  use  of  an  AC  supply  voltage  to  energise  the 
solenoid.  The  ADC  conversion  invariably  starts  before  the 
time  of  release  and  hence  one  could,  from  a  visual  inspection 
of  the  data  samples,  estimate  the  true  start  time,  to  the 
nearest  sampling  instant.  However,  the  error  in  this 
estimation,  which  depends  on  the  sampling  interval,  can 
significantly  influence  the  values  of  the  estimates  of  y*  and 
S*  . 


This  difficulty  was  overcome  here  by  fitting  the  exact 
analytical  solution  for  equation  (4)  to  the  experimental  data, 
using  an  optimisation  technique.  The  approach  is  illustrated 
by  Fig. 2.  If  tg  is  the  guessed  start  time,  and  tj  is  the  true 
start  time  then  there  are  two  cases:  (a)  tg  <  tj  (shown  in 
Fig. 2(a))  and  (b)  tQ  >  t?  (shown  in  Fig.2(b)).  In  both  cases, 
in  addition  to  the  unknown  parameters  y *  and  S*,  there  are  the 
unknown  variables  St  =  tj  tQ  and  A,  the  true  start 

amplitude.  It  is  therefore  required  to  minimise  the 

difference  between  the  analytical  solution  of  equation  (4)  and 
the  set  of  experimental  values,  Xj ,  I  =  0,1,..., n-1 ,  by 
optimising  the  values  of  St,  A,  y*  and  S*.This  is  achieved 
through  a  least-squares  fit  of  the  analytical  solution  to  the 
experimental  data. 

Thus,  if 

tA  =  iAt  -  St  (i  =  0,1 . n-1)  (7) 

then  the  analytical  solution  is 

X*  =  A  for  ti  <  0  (8) 

and,  for  damping  less  than  critical, 

*  -  y*xt.  y*x t. 

XA  =  Ae  [ sin ( — 2"^)  +  x  cos  ( — g-^)) 

for  tj  >  0  (9) 


where 


Similar  expressions  for  X<*  can  be  written  for  damping  equal 
to,  or  greater  than,  critical. 

A  least-squares  fit  is  achieved  by  minimising  the 
objective  function 

n=i 

U  =  jZq  (Xi  -  Xi*)2  (11) 


where  the  ranges  of  the  parameters  y*  and  S*  are  unrestricted. 
The  particular  optimisation  algorithm  adopted  here  is  a 
variant  of  the  PATSH,  a  direct  search  strategy  developed  by 
Whitney,  and  presented  by  Johnson  [22].  In  this  version  a 
progressively  centred  "shot-gun"  or  random  search  is  made  at 
the  end  of  each  direct  search  stage  in  an  attempt  to  avoid 
"hang-ups"  at  non-optimal  points.  In  the  present  application 
the  four-dimensional  optimisation  search  procedure  is  very 
much  simplified  due  to  the  absence  of  problem  constraints. 

Clearly  the  speed  of  convergence  of  the  optimisation 
method  depends  on  the  accuracy  of  the  initial  parameter 
estimates.  Here  the  parametric  identification  method, 

described  in  the  previous  section,  was  used  to  generate  the 
initial  estimates  of  y *  and  %*  ,  whilst  the  initial  estimates 
of  St  and  A  were  taken  to  be  0  and  1  respectively.  With  these 
start  values  the  optimisation  method  was  found  to  be  very 
rapid  and  efficient. 

5  EXPERIMENTAL  RESULTS 

For  each  test  condition  (i.e.  each  prescribed  value  of 
and  u0,  and  each  specified  direction  of  motion)  several 
(typically  five  to  seven)  individual  free-decay  records  were 
obtained.  Fig. 3(a)  shows  five  such  Individual  data  records 
for  a  typical  case  (eG  =  0.2,  wQ  =  80Hz,  radial  motion),  with 
a  sampling  interval  of  0.25ms.  Here  the  equi-spaced  samples 
have  been  linearly  Interpolated.  The  variations  in  the  data, 
within  this  set,  are  due  to  (a)  electronic  noise  in  the 
displacement  measuring  system,  (b)  analogue  to  digital 
quantisation  noise  (here  8  bit  resolution  was  employed)  and 
(c)  uncertainty  in  the  true  start  time,  as  discussed  earlier. 
However,  it  is  clear  that,  given  the  above  factors,  very  good 
repeatability  of  free-decay  data  was  obtained.  The  effect  of 
noise  could  be  substantially  reduced  by  averaging  the  set  of 
results  at  each  test  condition;  thus,  for  example,  the  effect 
of  averaging  the  data  shown  in  Fig.  3  (a)  is  shown  as  the  much 
smoother  decay  curve  in  Fig. 3(b). 

After  averaging,  at  each  test  condition,  the  data  were 
processed  using  the  techniques  described  briefly  in  section  4. 
Fig. 4  shows  some  typical  comparisons  between  the  fitted 
theoretical  curves  and  the  averaged  experimental  data,  for  the 
case  where  eQ  =  0,  wQ  =  180Hz  and  the  direction  of  motion  was 
radial.  For  data  acquired  at  the  relatively  large  sampling 
interval  of  0.5ms,  as  shown  in  Fig. 4(a),  it  is  evident  that 
the  parametric  Identification  technique  of  Ref. [21]  results  in 


a  rather  poor  fit  to  the  experimental  data  (note  that  here  the 
theoretical  curve  is  the  exact  solution  to  equation  (4), 
computed  using  the  estimated  values  of  y*  and  6* )  .  This 
discrepancy  arises  directly  from  the  use  of  a  relatively  large 
sampling  interval,  which  leads  to  (a)  significant  truncation 
errors  in  the  finite  difference  approximation  to  the  basic 
differential  algorithm  of  Ref. [21]  and  (b)  a  significant  error 
due  to  the  uncertainty  in  the  start  time.  However,  when  the 
optimisation  technique,  based  on  the  exact  analytical 
solution,  was  employed  as  a  second  stage  in  the  analysis,  then 
a  much  improved  fit  to  the  experimental  data  was  obtained,  as 
shown  in  Fig. 4(a).  It  can  be  expected,  therefore,  that  the 
parameter  estimates  obtained  from  the  optimisation  procedure 
are  significantly  more  accurate  than  the  original,  "identified 
parameters",  derived  by  the  method  of  Ref. [21]. 

As  the  sampling  interval  is  reduced  one  can  expect  that 
the  identification  method  will  improve  in  accuracy  and  that 
the  refinement  in  estimation  obtained  through  the  optimisation 
method  will  be  less  significant.  This  was  indeed  found  to  be 
the  case,  as  is  shown,  for  example,  by  the  results  in 

Fig. 4(b);  here,  for  the  same  test  conditions  as  in  Fig. 4(a),  a 
sampling  rate  of  0.25ms  was  used  to  acquire  the  data.  At  this 
higher  sampling  rate  it  is  evident  that  the  improvement  in 
fit,  through  the  second  level  of  analysis,  based  on 

optimisation,  is  much  less  significant. 

Table  2  shows  the  value  of  the  total  mean  square  error,  e, 
between  the  fitted  theoretical  decay  curves  and  the 

experimental  data  (for  the  optimisation  method  e  =  U/n,  where 
U  is  the  objective  function  defined  by  equation  (11));  this 
corresponds  to  the  results  shown  graphically  in  Fig. 4. 
Clearly,  for  At  =  0.5ms  the  value  of  e  has  been  reduced  by  up 
to  50%,  by  means  of  the  optimisation  technique,  whereas  at 
At  =  0.25ms  the  maximum  improvement  obtained,  over  the 

identification  technique,  is  about  12%. 

In  the  light  of  these,  and  other  similar  results,  it  was 
decided  to  process  all  the  averaged  decay  curves  using  the 
combined  identification-plus-optimisation  method  -  i.e.  using 
the  identification  procedure  to  provide  the  initial  guesses  of 
y*  and  S*  for  the  optimisation  calculation.  It  is  noted, 
however,  that  strictly  this  analysis  method  should  be  applied 
to  the  individual  free-decay  data,  prior  to  averaging,  and 
that  the  estimated  parameters  from  each  individual  data  record 
should  subsequently  be  averaged  to  produce  final  estimates  of 
y*  and  6*.  However,  experience  showed  that  the  values  of  y* 
and  S*  derived  from  the  individual  records,  when  averaged, 
were  very  close  to  the  values  deduced  directly  from  the 
averaged  data.  Some  typical  estimates  for  y*  and  S*  are  shown 
in  Table  3,  for  the  case  of  eQ  =  0.2,  wQ  *  80Hz,  and  radial 
motion  (see  Fig. 3  for  the  corresponding  experimental  data). 
It  is  seen  from  the  results  in  this  table  that  the  estimates 
derived  from  individual  tests  (see  Fig. 3(a))  show  a  low 
variability  and  that  their  averaged  values,  y*  =  0.553  and 
S*  =  0.898,  are  very  close  to  the  corresponding  values, 
y*  =  0.549  and  S*  =  0.907,  derived  by  analysing  the  averaged 
data  (see  Fig.3(b)).  All  the  results  presented  subsequently 
here  were  accordingly  derived  by  processing  the  averaged  data, 
at  each  test  condition,  since  this  involved  much  less 


computation  effort  than  the  processing  of  individual  data 
records . 

Figs. 5 (a)  to  (d)  show  some  typical  comparisons  between 
averaged  experimental  data  and  fitted  theoretical  curves, 
derived  by  the  identification-plus-optimisation  method.  These 
results  cover  the  frequency  range  50Hz  <  uD  <  180Hz.  In  all 
these  cases  a  very  good  degree  of  fit  is  achieved. 

In  Fig. 6,  the  final  estimates  of  damping  and  inertial 
coefficients  are  presented,  and  compared  with  the 
corresponding  short-bearing  theoretical  results,  computed  from 
the  expressions  given  in  Table  1.  The  damping  coefficients 
shown  here  are  the  non-dimensional  coefficients  per  land, 
defined  by 

Brr  =  brr  - ~  (12) 

2fmR$ 

and  similarly  for  Bss .  In  contrast  the  inertial  coefficients, 
crr  and  css  are  presented  in  dimensional  form,  and  have  the 
significance  of  hydrodynamic  mass,  mn  (in  kg),  for  both  lands. 
Results  are  presented  for  frequencies  of  vibration  in  the 
range  50  Hz  to  180  Hz.  It  is  noted  that,  in  the  experimental 
tests,  when  the  frequency  was  varied  all  the  other  variables 
were  maintained  the  same. 

The  experimentally  determined  coefficients,  Brr  and  Bss , 
in  Figs.  6  (a)  and  (b)  respectively,  are  seen  to  be 
significantly  greater  in  value  than  the  theoretical  values,  by 
a  factor  of  about  6,  at  low  eccentricities.  In  both  the 
radial  and  transverse  directions  the  coefficients  appear  to  be 
virtually  independent  of  frequency.  The  tendency  for  the 
experimental  coefficients  to  increase  with  increasing  eD  is 
reflected  in  the  theoretical  curves.  Moreover,  the 
experimental  values  of  Brr  and  Bss  are  very  similar,  at  low 
eccentricities,  which  again  is  in  line  with  the  theoretical 
predictions.  Thus,  qualitatively  there  is  general  agreement 
with  the  trends  predicted  by  the  short-bearing  theory;  this 
must,  however,  be  set  against  the  very  poor  quantitative 
agreement . 

It  is  noted  that,  in  previous  work  [15],  a  somewhat  closer 
agreement  between  measured  and  theoretical  values  of  Brr  and 
Bss  was  obtained,  although  experimental  values  were,  as  here, 
significantly  higher  than  the  theoretical  values.  In  Ref. [15] 
the  circumferential  groove  used  to  supply  lubricant  to  the 
squeeze-film  was  much  deeper  than  that  in  the  rig  used  in  the 
present  investigation.  It  is  possible,  therefore,  that  the 
geometry  of  the  supply  groove  is  an  important  factor  in 
determining  the  magnitude  of  the  damping  coefficients.  This 
is  a  factor  which  is  not  taken  into  account  in  the 
conventional  theory,  and  requires  further  investigation. 

The  experimental  inertial  coefficients,  shown  in  Fig. 6(c) 
and  (d) ,  for  the  radial  and  transverse  directions, 
respectively,  show  considerably  more  scatter  than  the 
corresponding  damping  results.  This  scatter  is  primarily  due 
to  the  fact  that  the  experimental  value  of  the  hydrodynamic 
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mass  (of  order  1  kg)  is  considerably  less  than  the  magnitude 
of  the  mechanical  first  mode  mass,  m  (of  order  10  kg).  Thus, 
in  computing  crr  from  equation  (6)  (and  similarly  css )  6*  is 
usually  close  to  unity  in  value.  This  results  in  errors  in 
the  estimation  of  s*  being  magnified,  by  a  factor  of  order  10, 
when  crr  and  css  are  computed. 

Bearing  in  mind  the  relatively  high  level  of  uncertainty 
in  the  measured  experimental  inertial  coefficients,  it  is 
possible  to  conclude  that,  with  the  possible  exception  of  the 
results  at  180Hz,  there  is  little  evidence  of  frequency 
dependency  in  the  coefficients.  The  most  significant  feature 
of  these  results  is  that  the  experimental  values  of  0^  and 
css  are  consistently  much  higher  than  the  corresponding 
theoretical  values,  the  discrepancy  being  typically  about  a 
factor  of  10.  In  this  respect  the  present  results  are  fully 
in  agreement  with  the  results  obtained  in  the  earlier 
investigation  of  Ref. [15].  There,  with  a  different  geometry 
of  squeeze-film  damper,  the  experimental  values  of  crr  and  css 
were  found  to  be,  typically,  about  an  order  of  magnitude 
higher  than  the  corresponding  short-bearing  theoretical 
values . 

In  assessing  these  results  it  should  also  be  borne  in  mind 
that,  as  mentioned  earlier,  each  value  of  static  eccentricity 
was  achieved  in  two  different  ways;  i.e.  with  the  Journal 
centre  above,  and  below,  the  bearing  centre.  In  the  radial 
direction,  the  direction  of  pull,  prior  to  release,  was  always 
upwards;  thus  differences  between  results  obtained  in  these 
two,  nominally  similar,  positions  can  be  attributed  to 
non-linear  effects,  arising  from  the  use  of  a  small,  but 
finite,  initial  start  amplitude,  d*  (typically  d*  «  0.2).  The 
maximum  discrepancies  encountered  were  for  the  case  of  uD  = 
50Hz,  and  high  values  of  c0 ;  here  the  two  values  of  crr  were 
found  to  differ,  by  up  to  100%,  for  eb  =  0.3  to  0.5;  these 
particular  results  must  therefore  be  regarded  as  unreliable. 

For  the  main  body  of  the  experimental  results  the 
vibration  amplitudes  were  well  within  the  linear  regime  of  the 
squeeze-film.  This  was  confirmed  by  performing  tests  at 
various  start  amplitudes.  When  normalised  to  unit  initial 
displacement,  the  free  decay  curves  were  found  to  collapse, 
for  different  values  of  d* ,  up  to  0.2,  indicating  that  a  small 
motion  (i.e.  linear)  model  should  be  applicable.  This  is  also 
confirmed  by  some  results  presented  in  Ref. 15,  where  the 
linear  regime  of  a  squeeze-film  was  explored  in  more  detail. 
As  one  would  expect  from  physical  reasoning,  the  degree  of 
collapse  tends  to  deteriorate  as  eQ  Increases. 

Returning  to  the  overall  results  shown  in  Fig. 6,  it  can  be 
concluded  that  there  are  serious  inadequacies  in  the 
conventional  short-bearing  theory  for  computing  both  damping 
and  inertial  coefficients.  In  this  respect  it  is  worth  noting 
here  that  it  has  recently  been  shown  that,  due  to  inertial 
effects,  large  pressure  variations  can  occur  at  the  inlet  to  a 
squeeze-film  bearing.  This  "Bernoulli  effect",  which  is 
critically  dependent  on  the  details  of  the  inlet  groove 
geometry,  is  not  taken  into  account  in  the  conventional 
theory.  In  the  present  experiments,  where  a  low  viscosity 
lubricant  was  used,  and  no  end  seals  were  attached  to  the 


bearing,  the  lubricant  flow  through  rate  was  relatively  high. 
In  these  circumstances  one  can  anticipate  that  inlet  inertial 
effects  are  likely  to  be  significant,  and  may  well  account  for 
the  large  discrepancies  observed  here  between  the  measured 
force  coefficients  and  the  corresponding  conventional 
theoretical  estimates. 

6  CONCLUSIONS 

The  main  conclusions  may  be  summarised  as  follows: 

1.  Damping  and  inertial  coefficients  can  be  estimated 
accurately,  and  efficiently,  from  digitised  free-decay  data. 
A  combination  of  parametric  identification  and  optimisation 
was  found  to  result  in  theoretical  curves  which,  in  general, 
very  closely  matched  the  experimental  data. 

2.  Both  the  measured  damping  and  inertial  coefficients 
were  found  to  be  significantly  higher  in  magnitude  than  the 
corresponding  values  deduced  from  conventional  short-bearing 
theory.  In  the  case  of  inertial  coefficients  the  observed 
discrepancy  between  theory  and  experiment  supports  the  results 
of  an  earlier  Investigation  [8]. 

3.  Both  the  damping  and  inertial  coefficients  were  found 
to  be  virtually  independent  of  the  frequency  of  vibration, 
over  the  frequency  range  50Hz  to  180Hz. 
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NOTATION 


direct  fluid-film  stiffness  coefficients 
for  the  radial  and  the  transverse 
directions,  respectively 

cross  fluid-film  stiffness  coefficients  for 
the  radial  and  the  transverse  directions 

matrix  containing  the  fluid-film  stiffness 
coefficients  arr  .  .  .  ass 

direct  fluid  film  damping  coefficients,  for 
the  radial  and  transverse  directions, 
respectively 

cross  fluid-film  damping  coefficients,  for 
the  radial  and  transverse  directions. 

matrix  containing  the  fluid-film  damping 
coefficients  bss  .  .  .  brr 

radial  clearance  between  the  journal  and 
the  bearing 

direct  fluid-film  inertial  coefficients  for 
the  radial  and  transverse  directions, 
respectively 

cross  fluid-film  inertial  coefficients  for 
the  radial  and  transverse  directions 

non-dimensional  direct  fluid-film  inertial 
coefficients  defined  by  equation  (16) 

non-dimensional  direct  fluid-film  inertial 
coefficients  defined  by  equation  (22) 

non-dimensional  direct  fluid-film  inertial 

coefficients  defined  by  equation  (26) 

matrix  containing  the  fluid-film  Inertial 
coefficients  css  ....  crr 


D  diameter  of  the  journal 

e  eccentricity  of  the  journal 

fr 1 (e,z) ,fs’ (e,z)  functions  of  e  and  z,  defined  in  equation 
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non-dimensional  fluid-film  forces,  due  to 
inertial  effects,  in  the  radial  and 
transverse  directions,  respectively 

fluid-film  pressure,  arising  from  inertial 
effects 

non-dimensional  fluid-film  pressure, 
arising  from  inertial  effects 
( P '  =  p/pu2R2) 

journal  displacement  in  the  radial  and 
transverse  directions,  respectively 

axial  distance  within  the  film 

non-dimensional  axial  distance  (z  =  Z/L) 

acceleration  of  the  journal  in  the  radial 
and  transverse  directions,  respectively 

non-dimensional  acceleration  of  the  journal 
(ar '  *  ar/u2c,  os'  *  as/u2c) 

static  eccentricity  ratio  (e/c) 

attitude  angle 

absolute  viscosity  of  squeeze-film  fluid 

density  of  squeeze-film  fluid 

angular  position,  measured  from  position 
of  maximum  film  thickness 


frequency  of  vibration 


SUMMARY 


Inertial  coefficients  for  full  squeeze  film  bearings  are 
evaluated  theoretically  using  Smith's  differential  equation 
relating  fluid  film  pressure  to  journal  acceleration  [1].  The 
variations  of  the  non-dimensional ised  inertial  coefficients 
with  static  eccentricity  ratio  in  the  radial  and  transverse 
directions  are  compared  with  some  corresponding  values 
obtained  from  Reinhardt  and  Lund  [2]  and  Szeri  et  al  [3]. 

The  results  from  these  three  methods  show  good  agreement, 
especially  for  short  bearings,  (i.e.  bearings  with  low  values 
of  length-to-diameter  ratio).  However,  Smith's  approach  has 
the  advantage  of  computational  simplicity  and  leads  to  fairly 
simple,  asymptotic,  analytical  expressions  for  very  short,  and 
very  long,  bearings. 

1.  INTRODUCTION 

Squeeze-film  dampers  have  proved  a  particularly  useful 
device  for  controlling  the  stability,  and  synchronous 
response,  of  the  rotating  shafts  which  occur  in 
turbomachinery.  In  such  machinery  there  has  been  a  tendency, 
in  recent  years,  to  increase  size  and  rotational  speed,  and  to 
use  light  oils,  of  low  viscosity.  All  these  factors  have 
resulted  in  increases  in  the  value  of  the  gap  Reynolds  number 


Re  =  pu)c2  -  inertia  force 
r\  ~  fluid  viscous  force 

at  which  squeeze-film  dampers  operate  in  practice;  values  of 
Re  of  order  one,  or  greater,  are  in  fact  now  common.  In  this 
situation  one  can  expect  inertial  dynamic  effects  to  be 
significant.  Recent  theoretical  studies  [1-8]  and 
experimental  investigations  [9-12]  have  indeed  confirmed  the 
important  influence  of  fluid  inertia  in  squeeze-film  dynamic 
behaviour . 

For  small  motions  of  the  squeeze-film  journal,  within  its 
clearance  circle,  it  is  possible  to  linearise  the  equations  of 
motion  and  to  represent  the  dynamic  behaviour  of  the 
squeeze-film  in  terms  of  damping  and  inertial  coefficients. 
Theoretical  methods  of  evaluating  the  damping  coefficients,  by 
using  Reynolds  equation,  have  been  known  for  many  years  but  it 
is  only  recently  that  the  problem  of  evaluating  the  inertial 
coefficients  theoretically  has  received  much  attention. 

Two  theoretical  approaches  to  this  problem  have  been 
reported  fairly  recently  by  Reinhardt  and  Lund  [2]  and  Szeri 
et  al  [3].  In  ref  [2]  a  perturbation  technique  is  used  to 
predict  values  for  the  inertial  coefficients  of  a  bearing  with 
a  finite  length-to-diameter  (L/D)  ratio.  Numerical  values  of 
these  coefficients  are  obtained  for  various  values  of  static 
eccentricity  ratio,  eD ,  and  for  L/D  values  of  0.1,  0.5  and  1. 
In  the  treatment  of  Szeri  et  al  [3]  an  alternative  approximate 
analysis,  derived  from  the  full  Navier-Stokes  equations  for 
fluid  flow,  is  also  used  to  compute  values  of  the  inertial 
coefficients  for  different  values  of  eG  and  L/D.  Their 
analysis  requires,  at  one  stage,  the  assumption  that  L/D  is 


small . 


Both  the  theoretical  treatments  of  refs  [2]  and  [3]  are 
involved  and  lengthy;  the  resulting  computation  procedures  for 
evaluating  the  inertial  coefficients  are  correspondingly 
complex  and  time-consuming.  There  is,  however,  a  much 
simpler,  alternative  theoretical  approach,  first  proposed  by 
Smith  [1],  which  leads,  in  turn,  to  much  simpler  calculations. 
This  approach  results  in  a  fairly  simple  differential  equation 
for  the  pressure  field  in  a  lubricant  film,  as  the  result  of 
journal  acceleration;  in  form  this  equation  is  similar  to  the 
familiar  two-dimensional  Reynolds  equation  for  finite 
squeeze-film  bearings.  Smith  did  not  attempt  a  numerical 
solution  of  this  differential  equation,  with  a  view  to 
evaluating  inertial  coefficients  for  finite  L/D  values,  but 
for  the  asymptotic  cases  of  L/D  -*  0  (short  bearings)  and  L/D  -► 
•  (long  bearings),  he  was  able  to  deduce  simple  analytical 
expressions  for  these  coefficients.  However,  despite  it's 
importance,  Smith's  approach  appears  to  have  gone  largely 
unnoticed  by  workers  in  the  area  of  squeeze-film  behaviour. 

The  objective  of  this  paper  is  threefold:  firstly  to  draw 
attention  to  the  relevance  of  Smith's  work  in  the  context  of 
squeeze-film  dampers;  secondly  to  obtain  values  of  the 
inertial  coefficients,  for  finite  L/D  ratios,  by  numerically 
solving  Smith's  pressure  differential  equation;  and  thirdly, 
to  compare  these  numerical  values  with  the  corresponding 
values  given  in  refs  [2]  and  [3]. 

2.  THEORY 

2 . 1  Hydrodynamic  forces 

An  examination  of  the  fluid  reaction  forces  on  an 
oscillating  body,  moving  in  an  incompressible  fluid,  indicates 
that  they  depend  on  its  displacement,  velocity  and 
acceleration.  Assuming  that  these  forces  depend  only  on  the 
instantaneous  motion  in  two  orthogonal  directions  (here  taken 
to  be  the  radial  and  transverse  directions,  r  and  s 
respectively) ,  the  general  linearised  expresssions  for  the 


hydrodynamic  forces,  Fr 

and  Fs,  for  small 

movements  of 

body  are,  in  matrix  form, 

,  as  follows  [1] 

F  =  ax  +  bx 

+  cX 

(1) 

where  F  =  Fr 

X  =  xr 

(2) 

Fs 

*S 

a  =  arr 

ars 
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asr 
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b  =  br  j. 

brs 

(4) 

bSr 

bSs 

and  c  =  crr 

crs 

(5) 

csr 

css 

Under  full-film  conditions,  (i.e.  no  starvation,  or 
cavitation),  the  cross-terms  disappear;  in  other  words,  the 
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matrices  a,  b  and  c  are  diagonal.  Furthermore,  all  the  terms 
in  a  are  negligible  [1].  b  and  c  contain,  respectively,  the 
damping  and  inertial  coefficients. 

2 . 2  The  inertial  coefficients 


A  linearised  approach  to  the  evaluation  of  a  squeeze-film 
fluid  forces  allows  a  separate  evaluation  of  the  Inertia 
coefficients  which  arise  from  the  effect  of  the  journal 
acceleration. 

Smith  [1]  has  shown  that  if  a'r  and  a's  are  the 
non-dimensional  journal  accelerations  in  .the  radial  (a'r  = 
xr/u2c)  and  the  transverse  (a's  =  xs/u2c)  directions 
respectively,  the  non-dimensional  fluid  pressure  p'  arising 
from  inertial  effects,  is  governed  by  the  following 
differential  equation 


fet(1  +  eo  cos  e)  ]  +  (?)2(1  +  eocos  e) 


=  a  '  cos  e  +  a  '  sin  e  (6) 

r  s 


This  equation  is  similar  to  the  two-dimensional  Reynolds 
equation  for  the  pressure  distribution  due  to  viscous  effects 
alone.  Here,  e  is  the  angle  measured  from  the  attitude  angle 
(e  *  0  at  maximum  film  thickness),  z  is  the  non-dimensional 
axial  position  within  the  fluid-film,  and  eQ  is  the  static 
eccentricity  ratio  (refer  to  Fig.l). 

The  solution  to  equation  (6),  subject  to  appropriate 
boundary  conditions  is  a  pressure  field,  p(e,z)  which  depends 
linearly  on  ar'  and  as'  .  The  general  solution  is  of  the 
form 


p'(e,z)  =  f r ' ( e , z ) ar '  +  fs'(e,z)as'  (7) 

The  fluid  inertial  forces  F'r  and  F's  (defined  as  positive 
when  applied  as  compressive  force  between  the  journal  and  the 
bearing)  may  be  evaluated  by  suitable  integrations  of  the 
pressure  field.  Thus,  these  forces  are  given  by: 


-  -  j  j  p 1 ( e , z )  cos  e  de  dz  (8) 

F'  =  -  [  f  p 1 ( © , z )  sin  e  de  dz  (9) 

s  JAJ 


where  integration  is  over  the  whole  pressure  field. 


It  follows  from  the  form  of  equations  (7),  (8)  and  (9) 
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Fs  1  —  Cgr  1  otj- 1  +  Cgg  as '  (11) 

where  the  non-dimensional  inertial  coefficients  are  given 
by 
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rr 
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sr 
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cos  e  de  dz 

cos  e  de  dz 


sin  e  de  dz 


sin  e  de  dz 


(12) 

(13) 

(14) 

(15) 


2.2.1  Evaluation  of  inertial  coefficients  for  finite  bearings 

The  general  pressure  differential  equation  (equation  (6)), 
was  solved  numerically  using  a  five  point  finite  different 
approximation  for  the  differentials  of  the  equation  in  the  e 
(circumferential)  and  z  (axial)  directions  for  a  range  of 
length  to  diameter  ratios  and  static  eccentricity  ratios. 

The  large  set  of  coupled  linear  equations  in  nodal 
pressure  was  solved  using  an  iterative  technique.  Depending 
on  the  trial  length  to  diameter  ratio,  with  the  demarcation 
arbitrarily  set  at  L/D  =  3,  the  nodal  pressure  field  was  set 
initially  to  the  Smith  analytical  values  for  either  short  or 
long  bearings  (the  short  bearing  solution  is  given  here  as 
equation  (18)).  The  iteration  sequence  was  terminated  when  the 
percentage  change  in  all  nodal  pressures  during  one  iteration 
was  less  than  0.1. 

It  was  assumed  that  the  bearing  consisted  of  two  Identical 
lands,  of  length  L,  separated  by  a  central  circumferential 
groove,  through  which  the  lubricant  was  supplied.  It  was  also 
assumed  that  the  lubricant  was  free  to  discharge  to 
atmospheric  pressure  at  the  bearing  outlets.  The  fluid  film 
pressure  arising  from  inertial  effects  was  taken  to  be  zero  at 
both  the  supply  groove  and  at  the  discharge  edges. 

Using  the  non-dimensional,  nodal,  pressure  field  and  with 
integration  over  the  full  land  width  axially,  and  over  2v  (for 
a  full  film)  circumferentially,  equations  (8)  and  (9)  were 
integrated  numerically  using  a  variant  Simpson's  rule  to 
determine  the  non-dimensional  fluid-film  forces  due  to 
inertial  effects.  Equations  (10)  and  (11)  were  then  used  to 
determine  the  non-dimensional  direct  fluid-film  inertial 
coefficients,  crr '  and  css'.  For  all  the  investigated 
geometries,  repeated  numerical  trials  were  made  with  finer 
finite  difference  grids  until  stable  inertial  coefficient 
values  were  obtained.  Generally  more  grid  points  were  needed 
to  achieve  stable  coefficient  values  for  long  bearings. 

It  should  be  noted  that  crr '  and  css '  can  be  converted  to 
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dimensional  form  as  follows: 


R3L  ' 

crr~  ^  c  crr 


R3L 
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where  crr  and  css  have  the  physical  significance  of  equivalent 
hydrodynamic  masses.  For  comparison  with  the  theoretical 
results  of  other  workers  and  with  Smith's  short  and  long 
bearing  analytical  solutions,  the  inertial  coefficients  were 
further  normalized  according  to  equations  (21)  or  (26) 
depending  on  the  trial  geometry  length  to  diameter  ratio.  The 
demarcation  length  to  diameter  ratio  was  taken  arbitrarily  at 
3.  The  comparisons  with  the  results  of  Reinhardt  and  Lund, 
and  Szeri  et  al,  thus  all  fall  in  the  short  bearing  category. 

3.  COMPARISON  OF  RESULTS 

3 . 1  Short  bearings 

Returning  to  the  general  differential  pressure 
relationship  given  by  equation  (6),  it  is  observed  that,  for  a 
short  bearing  (i.e.  L/D  <<  1)  the  pressure  gradient  in  the  z 
direction  will  dominate  over  that  in  the  circumferential 
direction.  In  these  circumstances  equation  (6)  may  be 
approximated  by 


(R)2  (1  +  cQcos  e) 


a2p  ' 


cos  e  +  a  '  sin  e  (17) 


On  integration  this  equation  gives  the  result 
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The  coefficients  in  equations  (12)  and  (15)  may  now  be 
evaluated  analytically;  the  result  is 
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For  the  present  purposes  it  is  convenient  to  introduce  the 
non-dimensional  coefficients 
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Equations  (19)  and  (20)  may  then  be  rewritten  as 
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for  the  non-dimensional  coefficients  crr"  and  css" .  In  the 
limiting  case,  as  cQ  -*  0,  crr"  =  css"  =  1. 

The  numerically  evaluated  inertial  coefficients  crr  '  and 
css '  for  finite  bearings  (see  section  2)  are  conveniently 
converted  into  the  non-dimensional  coefficients  crr"  and  css”, 
respectively,  for  the  purposes  of  comparison  with  the 
short-bearing  results  given  by  equations  (22)  and  (23). 
Clearly,  one  can  expect  the  inertial  coefficients  evaluated 
from  the  finite  difference  solution  to  approach  the  values 
evaluated  from  equations  (22)  and  (23)  for  low  values  of  L/D 
ratio  in  both  the  radial  and  transverse  directions. 

In  Figs.  2  to  6  the  variation  of  crr"  and  css " , 
respectively,  with  eQ  are  shown,  for  L/D  values  ranging  from 

0.1  to  2.0.  In  each  case  the  result  obtained  by  a  numerical 

solution  of  equation  (6)  is  compared  with  the  asymptotic 

results,  for  L/D  -♦  0,  as  given  by  equations  (22)  and  (23). 

Comparisons  are  also  made  with  the  results  presented  by 
Reinhardt  and  Lund  [2]  (values  given  for  L/D  =  0.1,  0.5  and 

1.0)  and  Szeri  et  al  [3]  (values  given  for  L/D  *  0.1,  0.2, 

0.5,  1.0  and  2.0). 

It  Is  observed  that,  for  L/D  <=  0.1  and  0.5,  the  values  of 
the  coefficients  obtained  from  the  different  methods  compare 
reasonably  well,  especially  for  the  transverse  direction. 
However,  for  the  same  L/D  value,  the  coefficient  values 
obtained  from  ref. [2]  are  slightly  higher,  in  the  radial 
direction.  The  agreement  tends  to  deteriorate  as  L/D 
Increases,  as  the  results  for  L/D  =  1  demonstrate.  The 

results  of  Szeri  et  al  [3]  agree  very  closely  with  the  present 
results,  for  L/D  «=  0.1,  0.2  and  0.5;  figs.  2  to  4  show  that 

there  is  virtually  exact  agreement,  in  this  range. 

The  difference  between  the  present  results  and  Szeri  et 
al's  results,  tends  to  increase,  as  L/D  increases.  This  is 
seen  clearly  by  the  comparison  of  results  for  L/D  =  2.0 

(fig. 6).  This  discrepancy  is  not  surprising,  since  Szeri  et 
al  [3]  introduce  the  assumption  that  L/D  is  small,  at  one 

stage  in  their  analysis.  Their  results  for  L/D  >  1  are  thus 
likely  to  be  less  accurate  than  the  present  values. 

It  is  noted  that  the  coefficient  values  obtained  from 

refs.  [2]  and  [3]  were  found  by  reading  values  from  graphs; 
small  numerical  errors  are  inevitably  introduced,  in  this 
process . 

3 . 2  Long  bearings 

For  the  long  bearings  (L/D  >>  1),  the  pressure  gradient  in 
the  circumferential  direction  is  dominant  and  the  general 
pressure  differential  relationship  of  equation  (6)  can  be 

approximated  as 
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This  can  be  solved  analytically.  Following  the  procedure  in 
section  3,  one  obtains  analytical  expressions  for  the  inertial 
coefficients;  thus 
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Here  it  is  convenient  to  scale  the  coefficients  so  that  the 
theoretical  values  are  unity,  at  e0  =  0.  Accordingly,  the 
non-dimensional  coefficients 


*  _  c. 


are  introduced  here.  Thus,  from  equation  (25) 
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for  the  asymptotic  case  L/D  -*  ®.  For  large  L/D  values,  the 
present  numerical  values  of  the  inertial  coefficients  Crr  '  and 
css'  (see  section  2)  will  be  converted  to  crr*  and  css* ,  for 
comparison  with  equation  (27). 

Fig.  7  shows  the  variation  of  the  inertial  coefficients 
with  static  eccentricity  ratio  for  different  values  of  L/D 
ratios,  as  compared  to  the  approximate  long  bearing  solution 
for  L/D  -♦  ■»,  for  the  radial  and  transverse  directions. 

At  high  values  of  L/C  ratio  the  numerical  values  of  the 
inertial  coefficients  approach  the  long  bearing  approximation 
As  the  L/D  ratio  decreases,  the  difference  in  the  inertial 
values  becomes  greater  as  shown  by  fig. 7. 

4.  CONCLUSION 


It  has  been  shown  that  numerical  values  of  the  inertial 
coefficients,  for  both  the  radial  and  transverse  directions, 
can  be  obtained  in  a  simple  way  by  numerically  solving  a 
partial  differential  equation  for  the  pressure  field  generated 
by  journal  acceleration.  The  numerical  values  obtained  were 
found  to  converge  towards  the  results  obtained  from  simple 
asymptotic  expressions,  for  both  L/D  -♦  0  and  L/D  «.  The 
results  presented  here  enable  the  accuracy  and  range  of 
validity  of  these  asymptotic  results  to  be  assessed.  A 
comparison  with  some  results  obtained  by  the  more  elaborate 
analyses  of  Reinhardt  and  Lund  [2]  and  Szeri  et  al  [3] 
indicates  that  the  present  approach  yields  results  of 
comparable  accuracy,  but  with  much  less  computational  effort. 
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CONCLUDING  REMARKS 

This  investigation  has  revealed  that  the  damping  and 
inertial  coefficients  of  a  squeeze-film  bearing  are,  at  least 
in  some  circumstances,  substantially  higher  than  the  values 
predicted  from  conventional  hydrodynamic  theory.  Thus,  to 
enable  the  formulation  of  adequate  design  procedures  for  these 
devices,  it  is  essential  that  the  source  of  the  discrepancy 
between  theory  and  experiment  be  identified  and  rectified. 

To  resolve  this  problem  it  is  suggested  that  the  present 
investigation,  described  here,  be  extended  to  determine  those 
operating  parameters  which  have  a  major,  controlling  influence 
over  the  discrepancy  between  theory  and  experiment.  It  is 
likely  that  the  inlet  geometry,  and  the  viscosity  of  the 
lubricant,  play  an  important  role.  For  example,  in  the  work 
described  here  the  flow-through  rate  of  the  lubricant  was 
high,  due  to  the  use  of  a  low  viscosity  lubricant.  With  high 
axial  flow  velocities  the  inertial  effects  at  the  squeeze-film 
inlet  are  likely  to  be  important,  and  may  well  be  dominant  for 
a  short  bearing.  Clearly,  changing  the  inlet  goemetry,  or  the 
lubricant  viscosity,  would  have  a  pronounced  influence  of  the 
magnitude  of  the  inlet  dynamic  effect.  This  effect  is  not 
taken  into  account  in  present  design  methods. 


Such  an  investigation  would  be  most  effectively  carried 
out  using  a  forced  vibration  test  method,  in  conjunction  with 
parametric  identification  techniques.  From  forced  vibration 
test  results  one  could,  for  example,  discriminate  between  mass 
(or  inertial)  and  stiffness  coefficients,  and  could  identify 
damping  parameters  even  when  the  damping  factor  is  greater 
than  the  critical  value  -  i.e.  the  case  "overdamped" 
vibration.  Thus  one  would  obtain  results  for  a  wide  range  of 
lubricant  viscosities,  in  a  rapid  and  efficient  manner.  Such 
work  would  involve  the  development  of  new  software  for 
parametric  identification,  and  some  modifications  to  the 
experimental  rig.  For  example,  it  would  be  necessary  to 
incorporate  a  force  transducer,  and  to  simultaneously  capture 
two  channels  of  information  -  displacement  response  and  force 
input . 

In  parallel  with  this  work,  a  more  detailed  experimental 
investigation  of  the  lubricant  flow  within  a  squeeze-film 
would  be  highly  beneficial.  The  most  promising  avenue  for 
research  here  would  seem  to  be  the  measurement  of  dynamic 
pressure  variations.  Through  investigations  of  the  variation 
of  the  instantaneous  pressure  distribution  with  time  it  should 
be  possible  to  pin-point  the  source  of  discrepancy  between 
theory  and  experiment. 

Finally,  fundamental  theoretical  studies  of  the  flow 
within  squeeze-film  are  required,  which  take  proper  account  of 
inlet  and  outlet  effects. 
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